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ABSTRACT 

We report the serendipitous detection by GALEX of fast (<lday) rising ( ^ Imag) UV emission 
from two Type II plateau (II-P) supernovae (SNe) at z=0.185 and 0.324 discovered by the Supernova 
Legacy Survey. Optical photometry and VLT spectroscopy 2 weeks after the GALEX detections link 
the onset of UV emission to the time of shock breakout. Using radiation hydrodynamics and non- 
LTE radiative transfer simulations, and starting from a standard red supergiant (RSG; Type II-P SN 
progenitor) star evolved self-consistently from the main sequence to iron core collapse, we model the 
shock breakout phase and the 55 hr that follow. The small scale height of our RSG atmosphere model 
suggests that the breakout signature is a thermal soft X-ray burst (Apeak ~90A) with a duration of 
< 2000 s. Longer durations are possible but require either an extended and tenuous non-standard 
envelope, or an unusually dense RSG wind with M ^ 10"'^ yr"-"^. The GALEX observations miss 
the peak of the luminous (Mpuv ~ —20) UV burst but unambiguously capture the rise of the emission 
and a subsequent 2 day long plateau. The postbreakout, UV-bright plateau is a prediction of our model 
in which the shift of the peak of the spectral energy distribution (SED) from '--^lOO to ~1000A and 
the ejecta expansion both counteract the decrease in bolometric luminosity from ~10^^ to ~10^Lq 
over that period. Based on the observed detection efficiency of our study we make predictions for the 
breakout detection rate of the GALEX Time Domain Survey. 

Subject headings: supernovae: general — ultraviolet: ISM — radiative transfer — stars: atmospheres 
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1. INTRODUCTION 

In recent years, with the advent of sensitive 7-ray and 
X-ray detectors, high-energy radiation from stellar ob- 
jects has been observed routinely. This includes in par- 
ticular the long-soft 7-ray bursts (GRBs), associated at 
present with a subset of hydrogen- and helium-deficient 
(Type Ic) SNe (Woosley & Bloom 2006). Their dura- 
tion of > 2s suggests a compact progenitor with a ~ls 
light-crossing time, a Wolf-Rayet (WR) star. Observa- 
tions have revealed that there is in fact a continuum from 
7-ray to X-ray signatures in association with such SNe. 
Hence, some events are related to GRBs in association 
with relativistic ejecta, while others are associated with 
X-ray flashes and the SN shock breaking out into the pro- 
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genitor's dense wind, as in the Type Ic SN 2006aj and 
the Type lb SN 2008D (Campana et al. 2006; Waxman 
et al. 2007; Soderberg et al. 2008.; Li 2007, 2008). In 
such SNe, the subsequent thermal emission from the cool- 
ing stellar envelope has also been observed (e.g., SN Ic 
1998bw: Galamaetal. 1998; SN Ib/c 1999ex: Stritzinger 
et al. 2002; SNe 2006aj and SN 2008D). 

In contrast to Type Ib/c SNe, hydrogen-rich (Type II) 
SNe originate from main-sequence objects in the mass 
range 8-25 M0, with iron cores at the time of collapse 
buried under an envelope of ~6~2OM0. The envelope 
light crossing time is ^0.65 -R3 hr (i?3 is the progenitor 
radius in 10"^ R0) and varies from about 5 minutes for a 
blue supergiant (BSC) progenitor to 1 hr in the largest 
RSGs (Levesque et al. 2005). In Type II SNe, these 
two properties prevent a central-engine from powering 
relativistic ejecta and from producing their associated 
non-thermal 7-ray/X-ray radiation. Hence, a burst of 
high-energy radiation from a Type II SN is most likely 
associated with shock emergence at the surface of the 
progenitor and its immediate aftermath. 

The shock-breakout phenomenon has been studied the- 
oretically (Falk 1978; Klein & Chevaher 1978), with spe- 
cific attention to the BSC progenitor of SN 1987A (En- 
sman & Burrows 1992; Blinnikov et al. 2000) or Type 
Ib/c SNe (Blinnikov et al. 2002). When the radiation- 
dominated SN shock wave travels outward through the 
envelope and reaches regions of moderate optical depth 
(i.e., r of a few tens), a radiative precursor initiates the 
brightening of the object leading to a brief soft X-ray 
burst. This precursor starts when the radiative diffu- 
sion time (idiff ^ 3r_ffp/c; Hp is the atmospheric scale 
height) becomes shorter than the shock-travel time to 
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the surface, located at r ^1. The envelope structure and 
extent are two key ingredients controlling the properties 
of the radiative precursor. In the absence of a wind, the 
expected durations are ~ls in a WR star (Blinnikov et 
al. 2002), -100s in a BSG (Ensman & Burrows), and 
~1000s in a RSG (Klein & Chevaher 1978). A very 
optically-thick wind, with a scale height comparable to 
the stellar radius, may, however, dramatically increase 
the diffusion time, and consequently the precursor dura- 
tion. 

In this Letter, we report the serendipitous detection 
with GALEX of rising UV emission in two Type II-P 
SNe discovered by the Supernova Legacy Survey (SNLS), 
SNLS-04D2dc (also reported by Schawinski et al. 2008) 
and SNLS-06Dljd, and associate both events with the 
breakout phase. This contrasts with previous early UV 
observations of Type II SNe, i.e., the Type Il-pec SN 
1987A (Hamuy et al. 1988), the Type lib SN 1993 J 
(Schmidt et al. 1993; Richmond et al. 1994), and the 
Type n-P SN 2006bp (Immler et al. 2007; Dessart et 
al. 2008, hereafter D08), in which they were all fading 
in the UV, since the SN ejecta had already begun their 
long-term cooling phase through expansion and radia- 
tion. In §2, we describe the GALEX and SNLS obser- 
vations. In fJSl we present our numerical simulations of 
the breakout and the early post-breakout evolution in a 
standard RSG progenitor, and compare these results to 
observations. In 21 we conclude and make predictions 
for the detection rate by the GALEX Time Domain Sur- 
vey (TDS). Throughout this Letter we adopt Hq = 70 
km s"^ Mpc-i, r2M=0.3, and 17a=0.7. 

2. OBSERVATIONS 

We conducted a search for early UV emission in Type 
II SNe by matching SNLS SN candidates with UV 
sources in overlapping GALEX Deep Imaging Survey 
(DIS) (Martin et al. 2005) fields. We selected SN can- 
didates that had ~ 1.5 ks GALEX DIS observations 
within 30 days before the optical SN candidate discov- 
ery, yielding a sample of 293 SN candidates with 1 to 77 
GALEX observations within that time period. Two vari- 
able UV sources were detected with a 3" matching radius 
to the position of SN candidates SNLS-04D2dc (R.A. 
= lOh 00m 16.681s, decl. = +02°12'18.52" [J2000.0]) 
and SNLS-06Dljd (R.A. = 02h 27m 36.189s, decl. = 
-04°31'56.62" [J2000.0]). The GALEX photometry was 
measured using aperture magnitudes with rap — 6" 
and a centroid fixed to that of the highest signal-to- 
noise (S/N) detection, with an aperture correction to 
account for the fraction of energy enclosed in the aper- 
ture of rritot — iTT-ap = —0.15 and —0.23 mag in the 
FUV (Aoff=1539A) and NUV (Aoff=2316A), respectively 
(Morrissey et al. 2007). 

SNe SNLS-04D2dc and SNLS-06Dljd were discov- 
ered on 2004 March 16 and 2006 December 12 by 
the SNLS real-time difference imaging pipeline (Astier 
et al. 2006; Sullivan et al. 2006). Follow-up VLT- 
FORSl spectroscopy was obtained on 2004 March 19 (JD 
2,453,083.5) and 2006 December 22 (JD 2,454,091.5), 
corresponding to 21.2 and 23.4 days after the GALEX 
detection, respectively. The SNe were classified as Type 
II-P SNe based on the identification of broad Hi emission 
lines (most noticeably II/3 and II7 in SNLS-06Dljd, but 
also Ha in SNLS-04D2dc), and a —40 day plateau in op- 



tical brightness. The host galaxy redshifts are measured 
from nebular emission lines and the spectra and pho- 
tometry are fitted with templates of a late-type galaxy 
at z = 0.185 with m,/=19.6 for SNLS-04D2dc and a 
starburst galaxy at 2 = 0.324 with mi/=22.1 for SNLS- 
06Dljd. Both VLT spectra match well the observations 
of SN 1999em at ~17 days after breakout (Dessart & 
Hillier 2006), suggesting that the observed UV brighten- 
ing is contemporaneous to the breakout phase. 

The host galaxies are detected in the FUV and 
NUV with mFuv=23.42±0.07 (23.80±0.08) and mNuv= 
22.94±0.06 (23.51±0.07) for SNLS-04D2dc (SNLS- 
06Dljd) measured from deep co-added images con- 
structed from observations taken before the SN detec- 
tion with tcxp= 37.6 ks (64.0 ks). We show the light 
curves of SNe SNLS-04D2dc in Figure 1 and SNLS- 
06Dljd in Figure 2, including GALEX FUV and NUV 
photometry (after subtraction of the host galaxy con- 
tribution) and SNLS difference imaging photometry in 
the g' (Aeff=4872 A), r' (Aoff=6282 A), i' (Aofr=7776 A), 
and z' (Acff=8870 A) bands. Note the >12day gap be- 
tween the UV and optical observations of both SNe. 
In order to relate the UV and optical photometry, we 
create synthetic GALEX- and SNLS-filter light curves 
using the multi-epoch non-LTE synthetic spectra com- 
puted for the Type H-P SN 2006bp (d = 17.5 Mpc; D08) 
and adjust to the distance, redshift, and extinction of 
SNLS-04D2dc and SNLS-06Dljd. For each, we adopt 
the onset of UV emission as the breakout time, yielding 
JD 2,453,062. 25lo 25 (JD 2,454,068. 05lo 74) for SNLS- 
04D2dc (SNLS-06Dljd). For SN 2006bp,' we adopt JD 
2,453,834.2 (Quimby et al. 2007). The optical colors dur- 
ing the plateau of both SNLS SNe are well fitted with a 
reddening of E{B — V) —0.1. However, an offset of -1-0.3 
(-0.4) for SNLS-04D2dc (SNLS-06Dljd) was added to 
obtain the best match to the observations, suggestive of 
a difference in intrinsic brightness of these SNe compared 
to SN 2006bp or an uncertainty in the extinction. 

3. MODELS AND RESULTS 

We model the UV and optical light curves from 0.05 to 
55.6 hr after shock breakout in three consecutive steps. 
First, using the one-temperature Lagrangian radiation 
hydrodynamics code KEPLER (Weaver et al. 1978), the 
evolution at solar metallicity of a 15 main sequence 
star is computed at high resolution (AM ~ 10~^ Mq 
at the surface) until iron core collapse. At that time, 
the star is a RSG with a radius i?*=865Ro and a to- 
tal mass of 12.58 M0 (Woosley & Heger 2007). The 
star above a mass cut of 1.45 Mq is then artificially ex- 
ploded with a piston to yield an asymptotic SN ejecta 
kinetic energy of 1.2x10^^ ergs typical of a core-collapse 
SN.^^ We finally post-process a sequence of 20 KEPLER 
models with the non-LTE non-relativistic steady-state 
code CMFGEN (Hillier & Miller 1998; Dessart & HiUier 
2005a) . The solution for level populations and the radia- 
tion field at —10^ frequencies is, thus, iterated to conver- 
gence on the stationary hydrodynamical structure, but 

However, originally, about half of the energy deposited by the 
shock wave is in the form of heat and trapped-radiation (internal 
energy), so that at shock breakout, i.e., 1.47 days after the piston 
trigger, the bulk of the SN ejecta in our model has not yet reached 
its coasting velocity. 
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aUowing the temperature to relax in regions exterior to a 
Rosseland-mean optical depth of ~20. This ensures that 
our assumption of steady-state is consistent, i.e., the dif- 
fusion time at the RSG surface being much less than the 
SN age. This is the first time the phase contemporaneous 
to shock breakout is modeled with a non-LTE multi-angle 
multi-frequency radiative transfer code, which allows for 
both line and continuum sources of opacities and treats 
multiple species and ionization states, including Hi, Hei- 
II, Cii-v, Nii-v, Oii-Ovi, Neiv-viii, and Feiv~xi and 
corresponding to a total of ~1000 levels. 

The one-temperature treatment adopted in KEPLER 
forces the gas and the radiation to the same tempera- 
ture, and, thus, prevents our modeling of the radiative 
precursor. However, the small atmospheric scale height 
Hp ^0.01 Ri, in our RSG progenitor surface suggests 
that the precursor will last ~1000s, light-travel-time ef- 
fects smearing the observed signal over ~i?^,/c ^^20008. 
Longer durations are possible, but would require a larger 
atmospheric scale height, perhaps reaching Hp ^ R^,. 
This may be produce d by an extended low-d ensity en- 
velope, as adopted bv lSchawinski et all (|2008D . although 
this appears incompatible with the RSG envelope struc- 
ture computed with KEPLER, the 3D hydrodynamics 
simulations of RSG envelopes (Freytag et al. 2002), 
and the quantitative spectroscopic analyses of RSG 
stars (Josselin & Plez 2007). Alternatively, this may 
be produced by an RSG in a super- wind phase, i.e., 
an RSG with an unusually high mass-loss rate M ~ 
10~^ M0 yr~^ compar ed to the more standa rd values of 
a few xlO-^Mgyr-i ()Salasnich et al.lll999D . 

In Figure[3l3, we show in black the evolution of the tem- 
perature Tphot and radius i?phot, defined by ^ ^ KtdR= 
2/3, where Kt = Ka + a and Ka (cr) is the absorptive 
(scattering) opacity. Tphot reaches up to ~1.35xlO^K 
only ~0.5hr after breakout, associated with a brighten- 
ing phase at constant radius (light-travel time effects and 
the radiative precursor would yield a ~2000s signal in- 
stead). Note that at this time, the "RSG" has the surface 
temperature of a white dwarf and radiates at 10^^ times 
the rate of the Sun! The temperature, decreasing almost 
as quickly as it rose and reaching down to 40,000 K af- 
ter ~5 hr, dominates the bolometric variation since mass 
shells have hardly moved during this early phase. At 
55 hr after breakout, Tphot ~20,000K, and i?phot has in- 
creased by only a factor of 3.8 to ^2.3x 10^'* cm. KE- 
PLER predicts a steep density distribution at the pho- 
tosphere, scaling as l/i?^*"^ at breakout to l/i?'"^" 55 hr 
later, a property also inferred from line profile morphol- 
ogy in SN 2006bp by DOS. With electron scattering 
dominating the opacity, the emergent SN flux suffers di- 
lution (Dessart & Hillier 2005b), reflecting better the 
temperature at the thermalization radius _Rth (defined 
by Jj^^^ KiiKtdR= 2/3) rather than at the photosphere 
(Fig. 3 a). These temperatures differ by up to a factor 
of 2-3 at all times, and the SED is much harder than 
the blackbody distribution at Tphot j with a peak (blue 
curve in Fig. Hi) at -90 A (900 A) at 0.5 hr (55.6 hr) 
after breakout. The non-LTE synthetic spectra shown in 
Figure [2t illustrate these properties and reveal for the 
first time the line features predicted under these break- 
out conditions (these characteristics will be discussed in 
a future paper). 



In Figures [T] and [2l we show the model GALEX and 
SNLS light curves created from our synthetic spectra for 
SNLS-04D2dc and SNLS-06Dljd. The CMFGEN models 
(1.67 hr < t < 55.6 hr) are offset by —1.5 (—1.4) mag for 
SNLS-04D2dc (SNLS-06Dljd) to match the SN 2006bp 
template light curves {t >2 days) that were scaled to fit 
the optical observations, suggesting that our models un- 
derestimate the observed luminosity by a factor of ~4. 
We find no significant variability in the photon data be- 
tween the first and second halves of the first two 20 
minute observations during the rise of the UV emission 
from SNLS-04D2dc and SNLS-06Dljd, with Am< 0.5 
mag. Although the lack of significant intra-observation 
variability and the — 2 mag amplitude of variability pre- 
dicted by our models for the <lhr long breakout signa- 
ture implies that the GALEX observations do not co- 
incide with the peak of the UV burst, the onset of the 
fast brightening UV emission is suggestive evidence that 
we have detected the breakout phase and constrains the 
time of shock breakout to an accuracy of 0.35 (0.75) days 
for SNLS-04D2dc (SNLS-06Dljd), i.e. 7.1 hr (13.7 hr) in 
the SN rest frame. Within the errors, the subsequent 
2 day long UV plateau is well reproduced by the model, 
although we note the poor match to the noisy FUV data 
in the plateau of SNLS-04D2dc and to the FUV/NUV 
color of SNLS-06Dljd at —1 day. Future, better sampled, 
higher S/N UV light curves of breakout events in Type 
II-P SNe will provide more accurate constraints, foster 
a better understanding, and, thus, permit more definite 
conclusions on this ephemeral breakout phase and the 
RSG envelope structure. 

4. DISCUSSION AND CONCLUSION 

GALEX DIS serendipitously detected fast rising UV 
emission in two SNLS Type II-P SNe, which we asso- 
ciate with the breakout phase in a RSG star. The ob- 
served amplitude of the UV signal is smaller than that of 
the < 1 hr long breakout burst predicted by our models, 
which suggests that the observations did not capture the 
peak of the burst (we see no strong variability within the 
20 minute observations), although our theoretical peak 
fiux is likely overestimated by the one-temperature treat- 
ment in KEPLER. The important conclusion from these 
models is that the observed UV brightening by '^Imag 
is predicted to occur only during this early, breakout, 
phase and thus is a direct probe of the time of shock 
breakout in these SNe. Moreover, we both observe and 
predict that the immediate post-breakout in the UV is 
followed by a plateau phase of about 2 days, supporting 
a RSG- like progenitor (more compact progenitors, e.g., 
BSGs, cool much faster after breakout, as seen for SN 
1987A), and providing an interesting clock to gauge the 
time since breakout. In SN 2006bp, the SN was fading 
in the UV right from the start of data collection, and 
a time since explosion of 1.7 days was indeed inferred 
(DOB). The post-breakout UV plateau is the result of two 
competing effects, with the strong decrease in bolometric 
luminosity compensated by the shift of the SED to longer 
wavelength as the ejecta expand and cool (Fig. [3h). At 
maximum light, the peak of the SED in our models is at 
90 A (in the rest frame) , and so the UV rise reveals only 
the evolution of the tail of the SED, not that of the peak. 

The volumetri c rate of Type H-P S Ne is close to that 
of Type la SNe (|Dahlen et al. Il2004f) : thus, we use the 
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measured fraction of Type la SNe in the total sam- 
ple of SNLS SN candidates from Sullivan et al. (2006), 
fia ~ 0.16, to estimate the fraction of Type ITP SNe 
in our sample. If we assume that the shock breakout 
occurs 15 ± 3 days before the optical discovery of the 
SN, then the number of candidates in our study that 
had a probability of at least one GALEX observation 
during the 2 day UV plateau following shock breakout 
is nobs = 167. By contrast, the number of candidates 
that had a probability of having a GALEX observa- 
tion that coincided with a 1 hr long flash at the time 
of shock breakout is ^ 1. We now estimate the num- 
ber of Type II-P SNe observed during the UV-plateau 
phase as nn-p ~ fiariobs ~ 27, and get a detection ef- 
ficiency of fdet = 2/nii-p - 0.07. GALEX TDS wiU 
operate in parallel with the Pan-STARRS Medium Deep 
Survey (PS MDS), an optical time domain survey that 
will monitor 50 deg^ night"^ with a depth comparable 
to SNLS. GALEX TDS plans to monitor -7deg2 of the 
PS MDS with cadences ranging from 98.6 minutes (two 
consecutive orbits) to several days for ^^8 months of the 
year to a depth of mum ^23 mag. PS MDS is expected 
to have a SN la detection rate of rja ~ 23 month" ^ 
per 7deg^ field (J. Tonry 2008, private communication), 
which for our observed detection efficiency {fdet ~ 0.07) 
and rimonths = 8 implies a favorable detection rate of 



UV emission from the shock breakout phase of Type II- 
P SNe by GALEX TDS of r^et « rianmonthsfdet ~ 10 
yr ^ 
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Fig. 1. — GALEX and SNLS (inset) light curves {colored symbols) of SNLS-04D2dc (z =0.185) shown relative to the estimated breakout 
time [horizontal bar, labeled ise) and overplotted with synthetic magnitudes {colored lines; II-P Template: scaled synthetic light curves 
of SN 2006bp [D08]; SB Model: shock-breakout model). Horizontal error bars show the time coverage of photometry measured from 
coadded observations. Solid arrows show 95% upper limits. The GALEX FUV detector was temporarily not operational during the NUV 
observations around t ~ 3.5 days. 
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JD-2454068.05 

Fig. 2. — GALEX and SNLS (inset) light curves (colored symbols) of SNLS-06Dljd (z =0.324), shown relative to the estimated breakout 
time (horizontal bar, labeled Ssb): and overplotted with synthetic magnitudes (colored lines; II-P Template: scaled synthetic light curves 
of SN 2006bp [DOS]; SB Model: shock-breakout model). Solid arrows show 95% upper limits. 
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Fig. 3. — Left: KEPLER/CMFGEN predictions for the post-breakout evolution of the wavelength of the SED peak (blue line; in units 
of 500 A), together with the gas temperature (Tphot; dashed line) and radius (Rphot; solid line) at the photosphere (black line) and at 
the thermalization depth (red line). Right: Multiepoch postbreakout synthetic non-LTE (thick line) and diluted blackbody (thin line; 
Tbb ~ ^th) spectra, scaled to the luminosity distance of SNLS-06Dljd. 



